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ABSTRACT  Random  amplified  polymorphic  DNA  (RAPD)  diagnostic  bands  are  one  tool  used  to 
differentiate  cryptic  mosquito  species  in  the  Anopheles  albitarsis  Complex.  Monophyly  of  four  species 
(A.  albitarsis  Lynch-Arribalzaga,  A.  albitarsis  B,  A.  deaneorum  Rosa-Freitas,  and  A.  marajoara  Galvao 
&  Damasceno)  currently  identified  with  the  RAPD  technique  was  assessed  using  sequences  of  the 
cytochrome  oxidase  I  (COI)  mitochondrial  DNA  (mtDNA)  gene.  Maximum  parsimony,  maximum 
likelihood,  and  Bayesian  analyses  support  monophyly  for  A.  albitarsis  s.s.,  A.  albitarsis  B,  and 
A.  deaneorum.  Anopheles  marajoara ,  as  identified  by  RAPD  banding  patterns,  was  either  polyphyletic 
or  paraphyletic  in  all  phylogenetic  analyses.  The  phylogenetic  pattern  and  within-species  genetic 
distances  observed  in  A.  marajoara  suggest  the  existence  of  a  previously  unidentified  species  (species 
E)  in  northern  Brazil  and  Venezuela.  Diagnostic  RAPD  bands  were  unable  to  distinguish  between 
A.  marajoara  and  species  E,  probably  because  of  the  low  number  of  correlated  bands  used  to  identify 
species  and  weaknesses  of  the  RAPD  technique,  in  particular,  violations  of  the  untested  assumption 
of  homology  of  comigrating  bands.  A.  marajoara  (even  without  species  E)  is  paraphyletic  with  respect 
to  A.  deaneorum;  if  A.  deaneorum  is  a  separate  species  from  A.  marajoara,  then  A.  marajoara  may  consist 
of  two  or  more  species  in  Amazonian  Brazil.  Based  on  mtDNA  COI  sequences,  there  are  at  least  four 
phylogenetic  species  within  the  Albitarsis  Complex:  A.  albitarsis  s.s.,  A.  albitarsis  B,  A.  marajoara,  and 
species  E;  the  species  status  of  A.  deaneorum  is  ambiguous. 
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Random  amplified  polymorphic  DNA  (RAPD)  analysis 
(Welsh  and  McClelland  1990,  Williams  et  al.  1990)  is 
a  polymerase  chain  reaction  (PCR) -based  technique 
that  uses  a  single,  short  oligonucleotide  primer  of 
arbitrary  sequence  to  produce  amplified  product. 
RAPD-PCR  allows  a  more  random  sample  of  the  ge¬ 
nome  (Black  1993)  with  rapid  production  of  large 
amounts  of  genetic  information  (Zhivotovsky  1999). 
RAPD-PCR  has  been  a  widely  used  tool  among  species 
in  the  genus  Anopheles  for  population  genetics  studies 
(Manguin  et  al.  1999,  Posso  et  al.  2003)  and  to  differ¬ 
entiate  cryptic  mosquito  species,  including  Anopheles 
gambiae  Giles  (Wilkerson  et  al.  1993),  A.  minimus 
Theobold  (Sucharit  and  Komalamisra  1997,  Kengne 
et  al.  2001),  and  A.  dims  Peyton  &  Harrison  (Manguin 
et  al.  2002). 

Nevertheless,  there  have  been  technical  (Riedy  et 
al.  1992,  Muralidharan  and  Wakeland  1993,  Schierwa- 
ter  and  Ender  1993)  and  theoretical  problems  noted 
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with  RAPDs  (Black  1993,  Clark  and  Lanigan  1993, 
Lynch  and  Milligan  1994).  A  common  assumption  is 
that  comigrating  bands  represent  homologous  loci. 
Comparisons  of  homologous  bands  for  inferring  phy¬ 
logenetic  relationships,  however,  may  be  misleading 
because  amplified  fragments  of  similar  size  from  dif¬ 
ferent  species  may  not  be  derived  from  the  same  loci 
and  not  all  RAPD  fragments  within  the  same  ampli¬ 
fication  pattern  are  independent  (van  de  Zande  and 
Bijlsma  1995). 

Anopheles  ( Nyssorhynchus )  albitarsis  Lynch-Arrib¬ 
alzaga  is  a  Neotropical  species  complex  that  includes 
at  least  four  cryptic  species  (Rosa-Freitas  et  al.  1990; 
Wilkerson  et  al.  1995a,  b),  only  some  of  which  are 
documented  malaria  vectors  (Wilkerson  et  al.  2005). 
Anopheles  deaneorum  Rosa-Freitas  is  the  only  species 
said  to  be  consistently  identified  using  morphological 
characters  (Rosa-Freitas  et  al.  1998).  Because  of  the 
lack  of  morphological  differentiation  among  the  cryp¬ 
tic  members  of  the  complex,  species  identification  in 
the  field  frequently  relies  on  the  presumed  geograph¬ 
ical  distribution  of  a  species  (Rosa-Freitas  et  al.  1998) . 
The  evidence  that  members  of  the  Albitarsis  Complex 
occur  sympatrically  in  some  locations  (Kreutzer  et  al. 
1976;  Steiner  et  al.  1982;  Narang  et  al.  1993;  Wilkerson 
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Table  1.  Summary  of  collection  localities  for  A.  albitarsis  s.l.  and  outgroups  A.  braziliensis  and  A.  darlingi 


Map  ref. 

State /Province 

Locality  (coordinates) 

RAPD-PCR  species0 

1 

Brazil 

Amazonas 

Manaus  (2°  53  S,  60°  15'  W) 

A.  marajoara  (CIO) 

2 

Bahia 

Itaquara  (13°  26  S,  39°  56'  W) 

A.  albitarsis  B  (Bl,  B4,  B5) 

3 

Ceara 

Fortaleza  (3°  43'  S,  38°  30'  W) 

A.  albitarsis  B  (B3) 

4 

Mato  Grosso 

Matupa/Peixoto  de  Azevedo  (10°  23'  S,  54°  54'  W) 

A.  marajoara  (C7,  C8) 

5 

Para 

Itaituba  (4°  15'  S,  55°  59'  W) 

A.  deaneorum  (D2) 

A.  marajoara  (C2,  C3,  C4,  C5,  C6) 

6 

Timboteua  (0°  56’  S,  47°  20'  W) 

A.  albitarsis  B  (B6) 

7 

Ilha  de  Maraj6fc  (0°  46'  S,  48°  31'  W) 

A.  marajoara  (Cl) 

8 

Rondonia 

Ariquemes  (9°  56'  S,  63°  04'  W) 

A.  marajoara  (C9) 

9 

Guajara-Mirim"  (10°  50'  S,  65°  20'  W) 

A.  deaneorum  (Dl,  D5,  D6) 

10 

Roraima 

Boa  Vista  (2°  49'  N,  60°  40'  W) 

A.  marajoara  (Cll,  C12,  C13) 

11 

Santa  Catarina 

Macaranduba  (26°  35'  S,  48°  58'  W) 

A.  albitarsis  s.s.  (A3) 

12 

Sao  Paulo 

6  km  SW  Registro  (24°  37’  S,  47°  53'  W) 

A.  albitarsis  s.s.  (A2) 

Amapa 

Macapa  (0°  2'  S,  51°  11' W) 

A.  albitarsis  B  (B2) 

A.  braziliensis 

Amapa 

Macapa  (0°  2'  S,  51°  11'  W) 

A.  darlingi 

13 

Argentina 

Buenos  Aires 

Baradero1'  (33°  48'  S,  59°  30'  W) 

A.  albitarsis  s.s.  (Al,  A4) 

14 

Corrientes 

Corrientes  (27°  28'  S,  59°  53'  W) 

A.  albitarsis  s.s.  (A6) 

15 

Misiones 

Posadas  (27°  23'  S,  55°  53'  W) 

A.  albitarsis  s.s.  (A5) 

16 

Venezuela 

Portuguesa 

Rio  Boeono  (8°  50'  N,  69°  59'  W) 

A.  marajoara  (C14) 

a  Individual  specimen  number  in  parentheses. 
h  Type  locality  of  A.  marajoara  (Galvao  and  Damasceno  1942). 

6  Type  locality  of  A.  deaneorum  (Rosa-Freitas  1989). 
dType  locality  of  A.  albitarsis  ( Lynch- Arribalzaga  1878). 

et  al.  1995a,  b)  further  complicates  determining  the 
roles  of  each  species  in  malaria  transmission. 

Although  there  has  been  recent  progress  on  the 
identification  of  diagnostic  primers  from  other  molec¬ 
ular  markers  (ITS2;  Li  and  Wilkerson  2005),  species 
identification  within  the  Albitarsis  Complex  currently 
relies  on  diagnostic  RAPD  bands  as  described  in  Wilk¬ 
erson  et  al.  (1995a,  b),  who  used  19  decamer  primers 
and  found  four  species  across  Argentina,  Brazil,  Par¬ 
aguay,  and  Venezuela:  A.  albitarsis  s.s.,  A.  marajoara 
Galvao  &  Damasceno,  A.  deaneorum,  and  a  fourth 
known  as  A.  albitarsis  B  that  has  yet  to  be  formally 
described.  These  four  species  are  predominantly  al- 
lopatric,  with  sympatry  recorded  at  nine  of  32  collec¬ 
tion  localities  (Wilkerson  et  al.  1995a) .  A  recent  study 
of  the  four  members  of  the  Albitarsis  Complex  based 
on  a  region  of  the  white  gene  that  includes  both  exon 
and  intron  sequence  determined  that  only  An.  mara¬ 
joara  retains  the  fourth  intron,  which  can  provide  a 
straightforward  way  of  distinguishing  this  species  from 
the  other  three  (Merritt  et  al.  2005). 

The  objective  of  this  study  was  to  assess  the  RAPD- 
PCR  methodology  of  Wilkerson  et  al.  (1995a,  b)  as  a 
species  identification  tool  in  A.  albitarsis  s.l.  by  com¬ 
paring  the  species  identities  resulting  from  RAPD  di¬ 
agnostic  banding  patterns  with  those  from  an  inde¬ 
pendent  marker,  the  mitochondrial  DNA  (mtDNA) 
cytochrome  oxidase  I  {COT}  gene.  COI  sequences 
have  been  shown  to  be  useful  in  resolving  evolution¬ 
ary  relationships  among  closely  related  species  groups 
for  a  wide  range  of  taxa  (Avise  1994,  2000),  including 
insects  (Brown  et  al.  1994,  Lunt  et  al.  1996)  and 
anopheline  mosquitoes  (Krzywinski  and  Besansky 
2003).  The  monophyly  of  the  four  putative  species 
found  by  Wilkerson  et  al.  (1995a,  b)  was  tested  to 


determine  whether  they  form  phylogenetic  species  as 
defined  by  Nixon  and  Wheeler  (1990). 


Materials  and  Methods 

Mosquito  Collections  and  Species  Identification. 
Collections  (Table  1)  were  made  between  1991  and 
1998  from  human  landing  catches,  surfaces  (blood-fed 
females  resting  indoors),  or  Shannon  traps  with  col¬ 
lectors  and  propane  lanterns  as  attractants.  Because  of 
the  polytypic  nature  of  the  complex,  specimens  from 
or  near  the  respective  type  localities  are  included: 
A.  albitarsus,  Baradero  (Buenos  Aires),  Argentina 
(Lynch- Arribalzaga  1878,  Rosa-Freitas  and  Deane 
1989);  A.  deaneorum,  Guajara-Mirim  (Rondonia) ,  Bra¬ 
zil  (Rosa-Freitas  1989);  and  A.  marajoara,  Ilha  de  Ma- 
rajo  (Para) ,  Brazil  (actual:  Cachoeira  do  Arari;  Galvao 
and  Damasceno  1942).  A.  albitarsis  B  has  yet  to  be 
formally  described  (Wilkerson  et  al.  1995a,  b).  All 
samples  sequenced  were  from  progeny  broods  as  de¬ 
scribed  in  Wilkerson  et  al.  (2005),  and  individuals 
from  each  family  are  retained  as  morphological  vouch¬ 
ers  in  the  Smithsonian  Institution,  National  Museum 
of  Natural  History  (NMNH).  The  lone  sample  from 
Portuguesa,  Venezuela,  collected  by  J.E.C.  was  a  wild- 
caught  female  adult.  DNA  from  this  female  is  main¬ 
tained  at  —  80°C  at  the  Conn  Laboratory  at  the  Wads¬ 
worth  Center. 

The  morphological  basis  for  the  definition  of  the 
subgenus  Nyssorhynchus  is  from  Peyton  et  al.  (1992). 
All  species  within  the  Albitarsis  Complex  were  iden¬ 
tified  morphologically  in  the  adult  stage  based  on 
Linthicum  (1988)  or  by  using  the  methods  described 
in  Wilkerson  et  al.  (1995b). 
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Total  DNA  was  extracted  from  individuals  using 
either  the  phenol-chloroform  method  described  in 
Wilkerson  et  al.  (1993)  or  the  Puregene  DNA  isolation 
kit  (Gentra  Systems,  Minneapolis,  MN).  For  each 
specimen,  species  identification  within  the  Albitarsis 
Complex  was  determined  using  the  RAPD-PCR 
method  of  Wilkerson  et  al.  (1995a,  b).  Individuals 
were  labeled  according  to  Wilkerson  et  al.  (1995b) 
with  species  A  as  A.  albitarsis  s.s.,  R  as  A.  albitarsis  R, 
C  as  A.  marajoara,  and  D  as  A.  deaneorum. 

Amplification  and  Sequencing  of  Mitochondrial 
DNA.  The  entire  COI  sequence  was  amplified  using 
two  sets  of  conserved  primer  pairs.  The  N-terminus  of 
the  gene  was  amplified  using  TY-J-1460  ( annealing  site 
located  in  the  tRNA-Tyrosine  coding  region;  5'-TA- 
CAATTTATCGCCTAAACTTCAGCC-3';  Simon  et 
al.  1994)  and  a  modification  of  UEA-6  (5'-TTAATTC- 
CTGTAGGNACAGCAATAATTAT-3';  Zhang  and 
Hewitt  1997).  The  COOH-terminus  of  the  COI  was 
amplified  using  modifications  of  Cl-J-2195  (5'-TGAT- 
TYTTTGGTCATCCN GAAGT-3' ;  Simon  et  al.  1994) 
and  TL2-N-3014  (5'-AATGCACTAATCTGCCATAT- 
TAG-3';  Simon  et  al.  1994).  The  PCR  reactions  were 
performed  in  25- /al  volumes  on  a  RoboCycler  Gradi¬ 
ent  96  Temperature  Cycler  (Stratagene,  La  Jolla,  CA) , 
each  containing  =1/100  the  DNA  of  a  mosquito,  20 
pM  each  primer,  and  one  Ready-To-Go  PCR  bead 
(GE  Healthcare,  Piscataway,  NJ).  An  initial  denatur- 
ation  of  3  min  at  95°C  was  followed  by  35  cycles  of 
amplification  (1  min  at  95°C,  1  min  at  53°C,  and  1  min 
at  72°C) .  The  products  were  purified  using  CentriSpin 
40  columns  (Princeton  Separations,  Adelphia,  NJ). 
The  fragments  were  sequenced  in  both  directions 
using  ARI  Prism  RigDye  version  1.0  Terminator  cycle 
sequencing  kit  (Applied  Riosystems,  Foster  City,  CA) 
and  run  on  the  ARI  Prism  377  automated  sequencer 
(Applied  Riosystems). 

Nucleotide  sequences  were  compiled,  edited,  and 
aligned  by  eye  with  the  aid  of  Sequencher  3.0  (Gene 
Codes  Corporation,  Ann  Arbor,  MI).  Aligned  se¬ 
quences  were  imported  into  PAUP*  4.0  (Swofford 
2003).  Sequences  were  aligned  with  the  published 
COI  DNA  and  protein  sequence  for  A.  g ambiae  (Gen- 
Rank  accession  no.  L20934;  Heard  et  al.  1993)  and 
translated  into  amino  acids  to  ensure  that  the  fragment 
sequenced  was  not  a  nuclear  pseudogene.  Sequences 
have  been  deposited  into  GenRank  (accession  nos. 
DQ076204  -DQ076227,  DQ076229  -DQ076234) . 

Phylogenetic  Analyses.  Five  to  14  individual  mos¬ 
quitoes  per  RAPD-PCR  species  were  sequenced 
(Table  1).  Two  individuals  of  each  of  A.  darlingi  Root 
(GenRank  accession  nos.  DQ076235  and  DQ076236) 
and  A.  braziliensis  Chagas  (GenRank  accession  nos. 
DQ076237  and  DQ076238)  were  included  as  out¬ 
groups.  According  to  Linthicum  (1988),  A.  braziliensis 
is  the  sister  taxon  to  A.  albitarsis  and  is  the  only  other 
member  of  the  Albitarsis  group  of  the  Argyritarsis 
section  of  Nyssorhynchus.  This  relationship  was  sup¬ 
ported  by  Sallum  et  al.  (2000)  based  on  a  cladistic 
analysis  of  163  morphological  characters.  A.  darlingi  is 
a  member  of  the  Argyritarsis  group  within  the  Argy¬ 


ritarsis  section  of  Nyssorhynchus  (Linthicum  1988) 
and  was  used  as  a  more  distant  outgroup. 

Roth  maximum  parsimony  (MP)  and  maximum  like¬ 
lihood  (ML)  analyses  were  implemented  in  PAUP* 
4.0.  Recause  of  the  number  of  ingroup  taxa,  optimal 
trees  were  generated  using  the  heuristic  search  option 
with  the  tree  bisection-reconnection  (TBR)  algo¬ 
rithm  chosen  for  branch-swapping  and  multiple  trees 
saved.  For  MP  analysis,  parsimony-uninformative 
characters  were  excluded,  and  100  random-taxon-ad¬ 
dition  replicate  analyses  were  carried  out  to  identify 
multiple  “islands  ’  (Maddison  1991)  of  most  parsimo¬ 
nious  trees  (MPTs).  For  ML  analysis,  an  appropriate 
substitution  model  and  model  parameter  values  were 
determined  using  the  computer  program  ModelTest 
3.0  (Posada  and  Crandall  1998).  This  program  gener¬ 
ates  and  statistically  evaluates  using  a  standard  likeli¬ 
hood  ratio  test  the  log  likelihood  scores  for  14  substi¬ 
tution  models.  Uncorrected  p  and  the  genetic  distance 
based  on  the  ML  model  were  calculated  using  PAUP* 
4.0. 

Nodal  support  was  determined  using  nonparametric 
bootstrap  proportion  (Felsenstein  1985)  and  Bremer 
support  (Bremer  1988).  Under  MP,  bootstrap  support 
values  were  generated  by  1000  pseudoreplicates,  with  10 
random-taxon-addition  replicates  per  pseudoreplicate, 
TBR  branch  swapping,  and  multiple  trees  saved.  Five 
hundred  pseudoreplicates  were  performed  under  ML, 
with  a  single  random-taxon-addition  starting  tree  per 
pseudoreplicate,  nearest-neighbor-interchange  (NNI) 
branch-swapping,  and  multiple  trees  saved;  NNI  branch¬ 
swapping  was  chosen  to  reduce  computational  time. 

Under  the  MP  framework,  the  Bremer  support 
value  or  decay  index  (Bremer  1988)  assesses  nodal 
support  by  measuring  the  tree  length  of  a  tree  con¬ 
strained  not  to  contain  a  particular  node,  subtracted 
from  the  unconstrained  most  parsimonious  tree 
length.  Index  values  measure  the  number  of  extra 
changes  required  on  the  tree  for  a  branch  to  be  lost; 
therefore,  higher  numbers  generally  indicate  greater 
support  at  a  node.  Decay  indices  were  calculated  for 
each  of  the  nodes  present  in  the  strict  consensus  of 
MPTs  using  the  AutoDecay  program,  version  4.0.1 
(Eriksson  1998). 

Mr  Bayes  version  3.0  (Huelsenbeck  and  Ronquist 
2000)  also  was  used  to  evaluate  optimal  trees  under  a 
Bayesian  phylogenetic  framework.  Program  default 
values  for  the  prior  probabilities  were  used.  The 
MCMC  algorithm  was  allowed  to  run  300,000  gener¬ 
ations,  and  sampled  every  100  generations  after  a 
burn-in  of  25,000  generations  (250  trees).  The  log 
likelihood  output  of  the  chain  was  plotted  to  confirm 
that  the  chain  had  stabilized  during  the  burn-in. 

Statistical  Tests.  The  optimal  trees  determined  using 
both  MP  and  ML  are  estimates  of  the  true  tree  and  may 
not  be  significantly  different  from  other  phylogenetic 
hypotheses.  Several  tests  under  both  MP  and  ML  were 
used  to  evaluate  whether  the  optimal  trees  produced 
by  COI  sequence  data  were  congruent  with  the  evi¬ 
dence  of  four  putative  species  identified  using  RAPD 
diagnostic  markers.  Constraint  trees  were  created 
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with  all  four  known  species  as  monophyletic  clades 
(within  species  branches  were  collapsed),  with  and 
without  sister  relationships  assumed  (A.  albitarsis 
s.s.lA.  albitarsis  B  and  A.  marajoaral A.  deaneorum ; 
Wilkerson  et  al,  1995a) . 

Under  the  MP  framework,  the  Templeton  test 
(Templeton  1983)  and  the  Kishino-Hasegawa  (KH) 
test  (Kishino  and  Hasegawa  1989)  were  used  to  de¬ 
termine  whether  MPTs  were  significantly  shorter  than 
the  constraint  trees.  The  Templeton  test  is  a  one-tailed 
Wilcoxon  signed-rank  test  that  compares  the  differ¬ 
ences  at  each  site  in  the  number  of  substitutions  re¬ 
quired  for  each  tree.  The  KH  test  is  a  f-test  that  com¬ 
pares  the  observed  differences  in  number  of  steps 
between  the  two  trees. 

Under  the  ML  criterion,  both  the  KH  and  the  Shi- 
modaira-Hasegawa  (SH)  test  (Shimodaira  and  Hase¬ 
gawa  1999)  were  used  to  compare  the  likelihood 
scores  of  the  ML  tree  with  suboptimal  trees  found 
under  constrained  searches.  Both  KH  and  SH  tests 
were  performed  using  resampling  estimated  log-like¬ 
lihood  bootstrap  (Goldman  et  al.  2000)  with  10,000 
replicates. 


Results 

Phylogenetic  Analyses.  Alignment  of  the  COI DNA 
sequences  produced  homologous  sequences  of  1,470 
bases  in  length  with  no  observed  indels.  Of  these,  297 
(20.2%)  sites  were  variable  and  243  (16.5%)  sites  were 
parsimony  informative.  The  strict  consensus  of  10 
most  parsimonious  trees  (MPTs;  tree  length,  528;  con¬ 
sistency  index  [Cl],  0.570;  retention  index  [RI], 0.793) 
is  shown  in  Fig.  1. 

Using  1,470  characters  from  the  ingroup,  the  best 
model  of  evolution  fitting  the  sequence  data  were 
HKY  +  I  +  T.  This  model  (Hasegawa  et  al.  1985) 
incorporated  unequal  base  frequencies,  a  transition: 
transversion  ratio  (6.9875),  proportion  of  invariant 
sites  (0.7691),  and  the  a  shape  parameter  of  the  T 
distribution  (0.9440).  Using  1,470  homologous  char¬ 
acters  for  both  the  ingroup  and  outgroups,  one  most 
probable  tree  was  produced  with  a  -Ln  likelihood 
score  of  5148.5398  (Fig.  2). 

For  Bayesian  analysis,  a  GTR  +  I  +  T  model  was 
used,  with  six  categories  of  rates  and  estimated  base 
frequencies.  The  strict  consensus  tree  is  shown  in 
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Fig.  1.  Strict  consensus  of  14  equally  parsimonious  trees  identified  by  COI  DNA  sequences.  Numbers  above  branches 
indicate  MP  bootstrap  proportions;  numbers  below  the  branches  indicate  Bremer  support  values.  Clades  I-IV  are  denoted 
in  italics.  Presumed  identities  (see  text)  are  as  follows:  A,  A.  albitarsis  s.s.:  B,  A.  albitarsis  B;  C,  A.  marajoara;  D,  A.  deaneorum ; 
and  E,  undescribed.  Outgroups  are  A.  braziliensis  and  A.  darlingi. 
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Fig.  2.  Single  tree  identified  by  maximum  likelihood  analyses  of  COI  DNA  sequences  under  the  HKY  +  I  +  T  model  of 
evolution.  Numbers  above  branches  indicate  ML  bootstrap  proportions.  Clades  I-IV  and  x,  y,  z  are  denoted  in  italics.  Presumed 
identities  (see  text)  are  as  designated  for  Fig.  1. 


Fig.  3.  Log  likelihood  scores  stabilized  (plot  not 
shown)  within  «=10,000  generations,  well  within  the 
burn-in  used  for  determination  of  posterior  probabil¬ 
ities. 

Monophyly  of  the  A.  albitarsis  s.s.,  A.  albitarsis  B, 
and  A.  deaneorum  species  clades  is  strongly  supported 
by  MP  (Fig.  1),  ML  (Fig.  2),  and  Bayesian  (Fig.  3) 
analyses.  Monophyly  for  A.  marajoara  as  identified  by 
BAPD  diagnostic  bands  is  not  supported;  this  species 
is  either  polyphyletic  or  paraphyletic  in  all  phyloge¬ 
netic  analysis. 

Among  all  haplotypes  within  the  Albitarsis  Com¬ 
plex,  mean  (range)  uncorrected  p  and  HKY  +  I  +  T 
( dHKY )  were  0.0360  (0.0170-0.0490)  and  0.0549 
(0.0200-0.0870),  respectively.  Mean  (range)  p  and 
dHKY  among  all  haplotypes  identified  as  A.  marajoara 
(C1-C14)  by  RAPDs  was  0.0335  (0.0007-0.0531)  and 
0.0517  (0.007-0.0978),  respectively.  Pairwise  genetic 
distances  within  and  among  each  putative  RAPD  spe¬ 
cies  and  A.  marajoara  clades  III  (C1-C10)  and  IV 
(C11-C14)  are  shown  in  Table  2. 

Topological  Constraint  Analyses.  Eight  optimal 
trees  (not  shown)  were  produced  under  the  con¬ 
straint  of  four  monophyletic  species  within  the  Albi¬ 


tarsis  Complex  (tree  length,  538;  Cl,  0.559;  and  RI, 
0.784).  When  hypothesized  sister  relationships  of 
A.  albitarsis  s.s. I  A.  albitarsis  sp.  B  and  A.  marajoara! 
A.  deaneorum  were  enforced  with  monophyletic  spe¬ 
cies  constrained,  34  optimal  trees  were  produced  (tree 
length,  597;  Cl,  0.601;  and  RI,  0.783).  The  10  uncon¬ 
strained  MP  trees  were  not  significantly  shorter  than 
either  the  eight  or  the  34  constrained  trees  under 
either  the  Templeton  or  KH  tests. 

Topological  constraints  under  the  ML  framework 
produced  the  same  tree  for  both  four  monophyletic 
species  with  and  without  sister  relationships  enforced 
(not  shown;  — Ln  likelihood,  5172.0699) .  This  tree  was 
significantly  less  likely  than  the  unconstrained  best 
tree  (difference  in  -Ln  likelihood,  23.5300;  PKH  test  = 
0.0168,  PSHtest  =  0.0127). 

Discussion 

COI  sequence  data  were  useful  in  resolving  some 
nodes  within  the  Albitarsis  Complex;  however,  deeper 
nodes  are  not  well  supported.  Among  MP,  ML,  and 
Bayesian  methods,  there  was  agreement  in  the  group¬ 
ing  of  individuals  into  at  least  four  distinct  monophy- 
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Fig.  3.  Strict  consensus  of  2,751  trees  identified  by  Bayesian  analysis  of  COI DNA  sequences  under  the  GTR  +  I  +  T  model 
of  evolution.  Numbers  above  branches  indicate  posterior  probabilities.  Clades  I-IV  and  x,  y,  z  are  denoted  in  italics.  Presumed 
identities  (see  text)  are  as  designated  for  Fig.  1. 


letic  lineages  with  one  paraphyletic  group.  Probably 
because  of  homoplasy,  relationships  between  haplo- 
types  and  lineages  lacked  resolution  under  MP,  and 
the  Templeton  and  KH  tests  were  not  able  to  detect 
significant  differences  between  the  MP  trees  and  con¬ 
straint  trees.  Relationships  among  the  lineages  were 
better  resolved  in  ML  and  Bayesian  analyses,  which 
incorporate  information  on  differing  substitution 
rates.  ML  demonstrated  significant  differences  be¬ 
tween  the  optimal  tree  and  the  constrained  tree  in 
statistical  tests,  with  the  optimal  tree  (showing  A. 
marajoara  as  polyphyletic)  being  more  likely  given 
the  COI  sequence  data. 

There  was  moderate  to  strong  support  (£83%)  for 
four  lineages  (clades  I-IV;  Figs.  1-3)  within  the  Albi- 
tarsis  Complex.  In  all  phylogenetic  analyses,  A.  albi- 
tarsis  s.s.  (A)  and  A.  albitarsis  B  grouped  together; 
however,  this  relationship  was  not  well  supported 
(bootstrap  or  posterior  probability  ranged  from  <50 
to  71).  Both  A.  albitarsis  s.s.  (A)  and  A.  albitarsis  B, 
identified  as  separate  species  by  fixed  RAPD  bands, 
form  strongly  supported  reciprocally  monophyletic 
clades  and  demonstrated  range  of  genetic  divergence 
(0.0259-0.0347  and  0.0336-0.0486  for  p  and  dHKY, 


respectively;  Table  2)  that  was  less  than  but  compa¬ 
rable  with  the  overall  range  observed  among  species 
(0.0170-0.0490  and  0.0200-0.0870  for  p  and  dHKY, 
respectively)  within  the  Albitarsis  Complex.  Other 
support  for  the  separate  species  status  of  A.  albitarsis 
B  can  be  based  on  a  chromosomal  species  concept 
(fixed  chromosomal  inversions;  Kreutzer  et  al.  1976, 
Steiner  et  al.  1982)  and  on  a  phenetic  species  concept 
(Steiner  et  al.  1982,  Rosa-Freitas  et  al.  1990). 

The  morphospecies  A.  marajoara  (based  on 
RAPDs)  seems  to  consist  of  at  least  two  phylogenetic 
species.  Either  polyphyly  or  paraphyly  of  the  group 
was  observed  in  all  phylogenetic  analyses.  In  statistical 
tests  under  ML,  the  optimal  tree  showing  the 
polyphyly  of  A.  marajoara  (C;  Fig.  2)  is  significantly 
more  likely  than  a  constrained  tree  with  A.  marajoara 
monophyletic.  Bremer  support  values  indicate  that 
eight  or  more  steps  would  be  required  for  A.  marajoara 
to  become  a  monophyletic  clade.  Moreover,  the  mean 
genetic  distance  within  A.  marajoara  (all  C  haplo- 
types)  (0.0335  and  0.0517  for  p  and  d^y,  respec¬ 
tively)  is  much  larger  than  the  within-species  genetic 
distance  values  (0.0074-0.0134  and  0.0079-0.0154  for 
p  and  dHKY,  respectively)  for  A.  albitarsis  s.s.,  A.  al- 


Table  2.  Average  (range)  pairwise  genetic  distance  within  and  among  CO  I  haplotypes  in  the  Albitarsis  Complex 
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bitarsis  B,  and  A.  deaneorum  (Table  2) .  There  are  four 
individuals  (C11-C14)  that  form  a  very  strongly  sup¬ 
ported  (>97%)  monophyletic  clade  (clade  IV;  Figs. 
1-3) .  The  mean  genetic  distances  (0.0467  and  0.079  for 
p  and  dHKY,  respectively)  between  A.  marajoara  of 
clade  III  (C1-C10)  and  clade  IV  are  comparable  with 
the  genetic  divergence  found  among  the  four  putative 
species  examined  in  this  study  (Table  2).  Given  the 
topology  of  all  optimal  trees  under  MP,  ML,  and 
Bayesian  frameworks,  the  strong  nodal  support,  and 
the  amount  of  genetic  divergence,  clade  IV  represents 
a  distinct  phylogenetic  species  within  the  Albitarsis 
Complex  (referred  to  as  species  E  for  the  remainder 
of  the  discussion). 

There  has  been  previous  genetic  evidence  in  sup¬ 
port  of  the  existence  of  a  different  form  (i.e.,  species 
E)  in  northern  Brazil  and  Venezuela  (Fig.  4) .  Based  on 
fixed  autosomal  and  X  chromosomal  inversions, 
Kreutzer  et  al.  (1976)  found  three  chromosomally 
differentiated  populations  within  A.  albitarsis  s.l.;  their 
form  C  found  in  Venezuela  and  Colombia  seems  to 
correspond  to  species  E  in  this  study.  In  addition, 
Rosa-Freitas  et  al.  (1990)  found  evidence  from  allo- 
zyme  loci  of  five  groups  based  on  Nei’s  genetic  dis¬ 
tance;  their  population  iii  found  in  Boa  Vista  ( Roraima 
state,  Brazil)  seems  to  correspond  to  species  E  in  this 
study.  Because  the  diagnostic  RAPD  banding  patterns 
for  A.  marajoara  and  species  E  are  indistinguishable, 
Wilkerson  et  al.  (1995a)  identified  all  samples  col¬ 
lected  from  Venezuela  as  A.  marajoara.  Despite  find¬ 
ing  some  morphological  differences  from  the  descrip¬ 
tion  ofLinthicum  (1988)  of  A.  marajoara,  Rubio-Palis 
et  al.  (2003)  have  only  identified  A.  marajoara  among 
all  A.  albitarsis  s.l.  specimens  collected  in  Venezuela 
because  they  used  RAPDs  for  species  identification. 

A.  deaneorum,  which  can  be  differentiated  from 
A.  albitarsis  s.s.  and  A.  marajoara  by  morphology  and 
fixed  RAPD  bands,  also  forms  a  well-supported  mono¬ 
phyletic  clade  within  the  A.  marajoara  lineage  and 
differs  genetically  from  the  other  three  putative  spe¬ 
cies  with  values  ranging  from  0.0170  to  0.0456  (p)  and 
0.0200-0.0808  (dHKY).  In  all  three  phylogenetic  anal¬ 
yses,  however,  A.  marajoara  (without  species  E)  is  not 
a  strongly  supported  reciprocally  monophyletic  clade 
separate  from  A.  deaneorum.  There  are  three  lineages 
within  the  A.  marajoara! A.  deaneorum  clade  (clade 
III)  under  both  ML  and  Bayesian  analyses  (clades  x, 
y,  z;  Figs.  2  and  3) ,  whereas  MP  shows  this  clade  to  be 
polytomous.  Anopheles  deaneorum  forms  a  strongly 
supported  monophyletic  clade  (clade  z)  under  ML 
and  Bayesian  analyses  and  a  moderately  supported 
monophyletic  clade  under  MP.  There  is  moderate 
bootstrap  support  under  ML  and  strong  support  in  the 
Bayesian  analysis  for  a  monophyletic  clade  containing 
Cl,  C3-5,  and  C7-9  (clade  x).  The  remaining  A.  ma¬ 
rajoara  haplotypes  (C2,  C6,  and  CIO)  form  a  strongly 
supported  monophyletic  clade  (clade  y)  under  MP, 
ML,  and  Bayesian.  The  relationships  among  clades  x, 
y,  and  z,  however,  are  unresolved  under  ML  and  enig¬ 
matic  in  the  Bayesian  analysis.  Mean  genetic  distances 
among  A.  marajoara  haplotypes  in  clade  III  and  be¬ 
tween  A.  marajoara  and  A.  deaneorum  are  congruent 
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Fig.  4.  Geographical  distribution  of  the  five  species  in  the  Albitarsis  Complex  determined  by  COl  DNA  sequences. 
Numbers  on  map  correspond  to  localities  in  Table  1. 


with  observed  genetic  distances  between  A.  albitarsis 
and  A.  albitarsis  B,  which  form  good  phylogenetic 
species. 

There  are  two  possible  explanations  for  the  ge¬ 
netic  patterns  and  lack  of  resolution  observed  in  the 
A.  marajoaral A.  deaneorum  clade.  The  first  is  that 
A.  marajoara  from  clade  III  may  consist  of  at  least  two 
phylogenetic  species.  The  monophyletic  clade  x 
(which  includes  Cl  from  the  type  locality  of  A.  ma¬ 
rajoara)  is  found  in  a  geographically  widespread  band 
across  Amazonia.  Some  A.  marajoara  haplotypes  from 
clade  y  (C2  and  C6)  are  found  in  sympatry  with 
haplotypes  from  clade  x  (C3-5)  in  Itaituba.  This  di¬ 
vergence  between  the  two  clades  may  have  occurred 
sympatrically  because  of  reproductive  isolation  or  it 
may  be  by  predominantly  allopatric  divergence  with 
secondary  contact  in  Itaituba. 

A  second  explanation  is  that  A.  deaneorum  and 
A.  marajoara  are  not  separate  species.  Interestingly, 
both  lineages  are  known  to  be  involved  in  malaria 
transmission  (Wilkerson  et  al.  2005),  and  our  findings 
lend  support  to  an  hypothesis  of  potential  introgres- 
sion  or  shared  ancestral  polymorphisms  for  genes  that 
would  confer  susceptibility  to  Plasmodi um.  A.  deaneo¬ 
rum,  as  described  by  Rosa-Freitas  (1989),  is  a  mor- 
phospecies  based  primarily  on  a  single  character: 
branched  outer  anterior  clypeal  hairs  in  the  larval 


stage.  This  character  seems  to  be  invariant  (R.C.W., 
unpublished  data),  and  separation  of  A.  deaneorum 
from  other  species  in  the  Albitarsis  Complex  is  further 
supported  by  allozyme  loci,  mtDNA,  and  hybridiza¬ 
tion  studies  (Rosa-Freitas  et  al.  1990,  Narang  et  al. 
1993,  Lima  et  al.  2004) . 

There  are  seven  diagnostic  RAPD  bands  that  are 
used  to  identify  A.  marajoara  (Wilkerson  et  al.  1995a, 
b);  however,  the  correlation  of  all  seven  bands  in  a 
single  sample  is  not  required  for  species  identification 
(Wilkerson  et  al.  1995b).  Commonly,  the  presence  of 
two  diagnostic  bands  is  used  to  identify  species.  Bands 
shared  by  A.  marajoara  and  species  E  may  be  because 
of  amplification  products  from  nonhomologous  loci  or 
amplification  from  areas  conserved  between  the  two 
species.  A  main  assumption  of  the  RAPD-PCR  tech¬ 
nique  is  that  similarity  of  band  size  is  a  dependable 
indicator  of  homology;  however,  nonhomology  of 
comigrating  RAPD  bands  has  been  documented  (Cog- 
nato  et  al.  1995,  van  de  Zande  and  Bijlsma  1995,  Riese- 
berg  1996) .  Among  220  pairs  of  comigrating  fragments 
in  wild  sunflowers,  Rieseberg  (1996)  found  only  79.1% 
useful  for  comparative  genetic  studies;  91%  of  the 
bands  exhibited  cross-hybridization  and/or  displayed 
congruent  restriction  fragment  profiles  suggestive  of 
homology;  however,  13%  of  these  loci  mapped  to  in- 
congruent  genomic  locations.  The  assumption  of  ho- 
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mology  of  similar  sized  fragments  for  RAPD-PCR,  in 
particular,  should  be  tested  to  determine  the  utility  of 
the  marker  in  phylogenetics  and  as  a  taxonomic  tool. 

Using  two  partial  sequences  of  mtDNA  COI  and 
ND4  genes  and  the  rDNA  ITS2  and  D2  expansion 
regions,  Wilkerson  et  al.  (2005)  used  total  evidence 
analysis  with  and  without  outgroups  to  examine  rela¬ 
tionships  among  members  of  the  Albitarsis  Complex. 
Their  study  and  the  present  one  both  detected  prob¬ 
lems,  using  one  or  more  analyses,  with  the  phyloge¬ 
netic  status  of  A.  marajoara  and  A.  deaneorum.  Despite 
the  strong  support  for  a  potential  fifth  species  in  our 
study,  we  recognize  that  our  findings  are  based  on  a 
single  gene  and  that  gene  trees  and  species  trees  do 
not  necessarily  detect  the  same  evolutionary  trajec¬ 
tory  (Presa  et  al.  2002).  To  better  understand  the 
species  status  of  A.  marajoara  and  A.  deaneorum,  we 
advocate  the  approach  in  Besansky  et  al.  (2003)  of  a 
locus-by-locus  assessment  of  members  of  species  com¬ 
plexes. 
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